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ABSTRACT 

We show how precision lensing measurements can be obtained through the lensing 
magnification effect in high redshift 21cm emission from galaxies. Normally, cosmic 
magnification measurements have been seriously complicated by galaxy clustering. 
With precise redshifts obtained from 21cm emission line wavelength, one can correlate 
galaxies at different source planes, or exclude close pairs to eliminate such contami- 
nations. 

We provide forecasts for future surveys, specifically the SKA and CLAR. SKA 
can achieve percent precision on the dark matter power spectrum and the galaxy 
dark matter cross correlation power spectrum, while CLAR can measure an accu- 
rate cross correlation power spectrum. The neutral hydrogen fraction was most likely 
significantly higher at high redshifts, which improves the number of observed galax- 
ies significantly, such that also CLAR can measure the dark matter lensing power 
spectrum. SKA can also allow precise measurement of lensing bispectrum. 

Key words: Cosmology: large-scale structure of Universe-radio lines: galaxies- 
galaxies: abundance 



1 INTRODUCTION 

Gravitational lensing measures the distortion of light by 
gravity originating from inhomogeneous distribution of mat- 
ter. The physics of weak gravitational lensing is clean, per- 
haps comparable to the primary CMB. To the first or- 
der approximation, it involves only general relativity and 
collision-less dark m atter dynamics. Ga s physics onl y en- 
ters at small scales dZhan fc Knoxll2004l ; IWhitdEool) and 
can be nulled by throwing awa y small scale lensing infor- 
mation jHuterer fe Whital2005l) . Though the prediction of 
weak lensing statistics is complicated by nonlinear evolu- 
tion of matter density fluctuation, high resolution, large box 
size N-body simulations are able to measure these statis- 
tics to high accuracy (e.g. lVale fc White! ll2uul) : lMerz et"ai~l 
(2005)). Thus, weak gravitational lensing is one of the most 
powerful and robust tools to constrain cosmology and study 
the large scale structure of the universe. 

Weak gravitat ional lensi n g has been detected through 
cosmic shear (see iRefregierl l j200 gj) for a recent review). 
Ongoing and upcoming large scale galaxy surveys such as 



CFHTLS 1 , DES 2 , LSST 3 , Pan-STARRS 4 and SNAP 5 will 
reduce the statistical errors of the cosmic shear measurement 
to < 1% level. At that stage, systematic errors such as point 
spread f unction, galaxy intrinsic alignment, seeing and ex- 
tinction iHoekstralEoollJarvis et al ]|2004t FVale et alJl20oJ : 
Ivan Waerbeke et al.l2004l) and errors in galaxy redshift mea- 
surement will be ultimate limiting factors. Given these pos- 
sible systematics, it is worthwhile to search for independent 
method to measure the gravitational lensing, for consistency 
check and, for likely better statistical and systematic errors. 
In this paper, we will show that the cosmic magnification of 
21cm emitting galaxies is a competitive candidate. 

Cosmic magnification is the lensing induced changes 
in galaxy number density. It introduces extra correla- 
tions in galaxy clustering jKaiserl Il992t IVillumsenl I199& 
iMoessner et alJll998tljain et alJl2003l) an d correlates galax- 
ies (q uasars) at widely separate redshifts llMoessner fc Jainl 
1998). Cosmic magnification contains as much informa- 
tion of cosmology and matter clustering as cosmic shear. 
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It has been robustly detected in quasar-galaxy lensing 
iScranton et al.l i2005t) and references therein). 

The measurement of cosmic magnification does not re- 
quire the accurate determination of galaxy shapes and is 
thus free of many systematics, such as point spread function 
and galaxy intrinsic alignment, entangled in the cosmic shear 
measurement. But it suffers from several obstacles. (1) It suf- 
fers stronger shot noise than cosmic shear measurement. For 
cosmic shear, shot noise comes from the intrinsic ellipticity 
of galaxies, which has dispersion (e 2 ) 1 / 2 ~ 0.3. The shot 
noise power spectrum is proportional to (e 2 )/N g ~ 0.1 /N g , 
where N g is the total number of galaxies. For cosmic mag- 
nification, shot noise comes from the Poisson fluctuation of 
galaxy counts and the shot noise power spectrum scales as 
1/Ng and is thus an order of magnitude larger than that 
in cosmic shear. (2) The signal of cosmic magnification is 
generally much smaller than intrinsic clustering of galaxies. 
Without precise redshift measurement of galaxies or quasars, 
the only method to remove intrinsic clustering of galaxies 
is to measure the cross correlation of foreground galaxies 
and background galaxies (or quasars). But its detection is 
severely limited by the number of high redshift galaxies, 
which are difficult to detect in optical band, and quasars, 
which are rare, comparing to galaxy abundance. Early mea- 
surements of cosmic magnification were often controversial. 
Even with relatively large sample of quasars and foreground 
galaxies from 2dF and SDSS, the measurement of cosmic 
magnification is still at its infancy and is only confirmed to 
^ 8cr confidence level. 

Besides the observational difficulties, the theoretical 
prediction of cross correlation strength is complicated by 
foreground galaxy bias, which is hard to predict from first 
principles. Without precise understanding of the galaxy bias, 
the power of cosmic magnification to constrain cosmology 
and matter distribution is severely limited, compared to cos- 
mic shear. 

Cosmic magnification measured from 21cm emitting 
galaxies is free of many obstacles entangled in the cosmic 
shear measurement and quasar-galaxy cosmic magnification 
measurement. As we will show in §5] upcoming radio surveys 
such as Square Kilometer Array (SKA) 6 can find ~ 10 8 -10 9 
Hi-rich galaxies in total and ~ 10 7 galaxies at z > 2 through 
the neutral hydrogen 21cm hyperfine transition line. With 
this large sample of galaxies, the measurement of cosmic 
magnification will enter the precision era, due to following 
reasons. (1) Radio observations are free of extinction. Since 
dust is associated with galaxies, extinction is correlated with 
foreground galaxies and thus is likely able to produce a false 
quasar-galaxy or galaxy-galaxy cross correlation signal. (2) 
Contrary to optical surveys, redshifts of these galaxies can 
be precisely determined by the redshifted wavelength of the 
21cm hyperfine transition line, with no added observational 
cost. Precision measurement of galaxy redshift is required for 
the precision prediction of lensing statistics. It also allows 
the luxury of removing close galaxy pairs, which is crucial to 
measure cosmic magnification in the presence of the strong 
auto correlation of galaxies. (3) At high redshift the magni- 
fication effect is enhanced, since galaxies at higher redshifts 
are more strongly lensed and have a steeper flux distribu- 
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tion. Combined with the large sample of galaxies, shot noise 
can be overcome. (4) The auto correlation function of cos- 
mic magnification can be measured to high accuracy. The 
prediction of the cosmic magnification auto correlation is 
free of galaxy bias prior and as robust as the prediction of 
cosmic shear power spectrum. 

The goal of this paper is modest, which is to demon- 
strate the feasibility to do precision lensing measurement us- 
ing 21cm emitting galaxies. Methods discussed in this paper 
are by no mean optimal and results presented are conserva- 
tive. Efforts toward optimizing lensing measuremen t is pre- 
sented in the companion paper | Zhang fc Penll2005l) . where 
we show that cosmic magnification of 21cm emitting galaxies 
can not only do better than of cosmic magnification optical 
galaxies, but also can likely do better than cosmic shear of 
optical galaxies. 

Throughout this paper, we adopt a flat ACDM uni- 
verse with Q m = 0.3, S1a = 1 — f2m, erg = 0.9, h — 0.7 
and in itial power index n = 1, as is consistent with WMAP 
result iSpergel et al.ll2003T) . We take Canadian Large Adap- 
tive Reflector (CLAR) 7 and SKA as our targets to forecast 
the ability of future radio 21cm surveys to measure cosmic 
magnification. The instrumental parameters and survey pat- 
terns of these two surveys have not been completely fixed 
yet. For CLAR, we adopt a system temperature T sys = 30 
K, effective collecting area A c g — 5 x 10 4 m 2 and field of 
view 1 deg 2 . For SKA, we adopt the same value of T Bys and 
field of view and adopt A e g — 6 x 10 5 m 2 . 



2 DETECTING GALAXIES IN RADIO 
SURVEYS 

Though the universe is highly ionized, there are still large 
amounts of neutral gas p resent in galaxies. The typical HI 
mass is around ~ 10 9 Mq iZwaan et alll99^) . HI rich galax- 
ies appear in the radio band by emitting the 21 cm hyperfine 
line resulting from the transition of atomic hydrogens from 
spin 1 ground state to spin zero ground state. The sponta- 
neous transition rate is A21 = 2.85 x 10 _15 s _1 . 

Each emission line has negligible width. But since neu- 
tral gas has rotational and thermal motions, the integrated 
intrinsic line width is not negligible. It turns out that the line 
width of 21 cm emission is mainly determined by the rota- 
tion of HI gas, which has a typical velocity ~ 100 km/s. The 
Doppler effect by thermal motions cause a velocity width 
w ~ Cy/k B T/m H < 10^T/10 4 K km/s. The actual emis- 
sion line width is determined by many parameters, such as 
the total mass of galaxies, redshift and inclination angle of 
the HI disk. For simplicity, we assume that the combined 
intrinsic line width is w = 100 km/s. The choice of w affects 
the prediction of detection efficiency of HI emitting galaxies. 
But the results shown in this paper should not be changed 
significantly by more realistic choice of w. 

Since the velocity (frequency) resolution of radio sur- 
veys can be much higher than the intrinsic line width, the ob- 
served flux could have a nontrivial dependence on the band- 
width, whose modeling requires detailed description of HI 
distribution in galaxies. To avoid this complexity, we adopt 

7 http://www.clar.ca/ 
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bandwidth to be larger or equal to redshifted 21cm line 
width, which changes from w at z to w/(l + z) at z — 0, or 
in frequency space, to Af = wvni/c[l + z), where V21 = 1.4 
Ghz is the 21cm frequency. The total 21 cm flux of HI rich 
galaxies is 

32^21 N m E 2 i ,w V21 



S-> 



.gi + .92 4nD 2 L (2 



-(■ 



= 0.023mJy 



c 1+z' 



(1) 



lOOkm/s 

1O 1O M v x(z) ' w(l + z) 



Here, g 1 = 1 and 92 = 3 are the degeneracies of atomic 
hydrogen spin and 1 ground state. E21 = hi>2\ is the energy 
of each 21cm photon. iVm = Mm/mm is the total number 
of neutral hydrogen of a galaxy with total hydrogen mass 
Mm. Dl(z) = x( 2: )(l + z ) is the luminosity distance and x 
is the comoving angular distance. Ho = 100ft km/s/Mpc is 
the Hubble constant at present time. 

Instrumental noise scales as Af^ 1 ^ 2 , where Av is the 
bandwidth. So a larger bandwidth helps to beat down in- 
strumental noise. But if the bandwidth is larger than the 
21 cm line width, the signal is diluted and scales as 1/Av. 
Thus, the highest signal to noise ratio is gained when the 
bandwidth Av is equal to the redshifted bandwidth of inte- 
grated 21cm emission line, wv2\jc(\ + z). The system noise 
per beam is 



\/2T s , 



0.032mJy r - 5Xl0V2 



(2) 
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Here, 
Vc = 



n c is the correlator efficiency, which is adopted as 
0.9. The instrumental beam area is At ~ \ 2 /A e g. 
For CLAR, A b ~ 10 2 (v 2 i / v) 2 ■ For SKA, A b ~ l' 2 (v 21 /v) 2 . 
21cm emitting regions have typical size ~ 30 kpc/h or ~ 1 
at z ~ 1. Thus, the size of 21cm emitting regions is much 
smaller than the beam size. The observed flux is thus the 
total flux of each galaxies. In this case, the calculation of 
the detection threshold is straightforward. 

If we choose those peaks with flux above nS sys (n-cr 
selection threshold), the minimum HI mass selected is 



n x 1.39 x 10 lu M o 



c/H Q 



(l + z 



3/2 



w hour T sys 5 x 10 4 m 2 
lOOkm/s t 30K 



(3) 



We assume that the HI ma ss function follow s the 
Schechter function found at z = dZwaan et al]ll997h 



n(M,z)dM = n (z) 



M 



M,{z) 



exp 



M 



M*{z) 



dM . (4) 



We fix 7 = 
0) = 10 9 ' 55 /i 



1.2. IZwaan et alJ 

~ 2 Mq and no(z — 



found that M*(z = 
0.014/i 3 Mpc -3 . There 



is little solid measurement of no (z) and M,(z) other than 
at local universe. But for this form of the mass function, 
there exists a tight relation between Ohi, the cosmological 
neutral hydrogen density with respect to the present day 
critical density, and noM*: 



Qmh = 2.1 x 10" 
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Figure 1. The predicted abundance of 21cm emitting galaxies 
above 4cr detection threshold of CLAR 5 year survey. The solid 
lines are the cumulative number distribution JV(< 2) and the 
dashed lines are the differential distribution dN/dz. The total 
number of observed galaxies is of the order 10 6 and the number 
of galaxies at z > 1 is of the order 10 5 . We assume no evolution 
in the mass function. If realistic evolution model is considered, 
the total number of galaxies at z > 1 can increase by a factor of 
5 or more. 



Observations of damped Lyman-a systems and Lyman- 
q limit systems measure Qui from z = to z ~ 4. Com- 
bining Eq. and these observations, one can put con- 
straints on no and M*. These observations found that 
!Qhi increases by a factor of 5 toward z ~ 3 and 



then decreases toward h i gher r edshift (Zwaaii^^l 



toward nig 

Rao fc Turnshe k| l| 2 000D ; jgtorrie-Lombardi fc Wolfe 



1997 
(2000) 



P^roux^tHuTT^OOsl) and data compilation in IPeToux^t^rt] 
l2^03T? T l^agamme et alJ <l2004f) i. Either the increase of no 
or M» increases the detectability of 21cm emitting galax- 
ies. Thus, estimations based on the assumption of no evo- 
lution should be regarded as conservative results. We will 
show that even in this conservative case, future radio sur- 
veys such as CLAR and SKA still allow precise measure- 
ment of galaxy-galaxy lensing (SjJJl and in the case of SKA, 
the precise measurement of lensing power spectrum (33]&|3J 
and lensing bispectrum ([0. If we adopt evolution models 
implied by and consistent with observations, even CLAR is 
able to measure the lensing power spectrum 

The number of galaxies detected depends on the se- 
lection threshold n. If we choose 4cr detection threshold 
(n = 4), CLAR could detect ~ 10 6 -10 7 galaxies in a 5-year 
survey. A deeper survey (smaller sky coverage) detects a 
larger fraction of high redshift galaxies. But since the survey 
volume is smaller, the total number of high redshift galaxies 
is not necessarily higher. A survey area around 100 deg 2 is 
optimal to detect high z galaxies. For a 160 deg 2 survey area, 
~ 10 6 galaxies at z > 1 can be detected. SKA is about 10 
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Figure 2. Similar as Fig.Q but for SKA 5 year survey. The total 
number of observed galaxies is of the order 10 8 and the number 
of galaxies at z > 1 is of the order 10 7 . We assume no evolution 
in the mass function. If a realistic evolution model is considered, 
the total number of galaxies at z > 1 can increase by a factor of 
5 or more. 



times more sensitive than CLAR and can detect two orders 
of magnitude more galaxies. For a 1600 deg 2 survey area, 
~ 10 s galaxies at z > 1 and ~ 10 7 galaxies at z > 2 can 
be detected. As a reminder, these estimations are extremely 
conservative. The number of galaxies detected at z > 1 can 
be enhanced by a factor of 5 or more if the evolution effect 
is considered. 

Some peaks above the selection threshold are caused by 
noise. N no i Be , the number of false peaks, has strong depen- 
dence on the detection threshold. If we choose n = 1( la 
detection), the number of false peaks of noise is 



A no 



Anf sky v Erfc( 1.0/^2) 10 i 
7. ^ 10 /sky 



'2 



.4 



pixel 



-4 6 



(6) 



For CLAR deep surveys which cover / s k y ~ 1% of the sky, 
A^noise < 10 7 , which is still less the number of detected 
galaxies above la threshold. For SKA with several thousand 
square degree sky coverage, AWse above la is ~ 10 8 , which 
is slightly less than the total number of galaxies above la 
threshold. In CMB measurements, signal-to-noise per pixel 
is often chosen to be 1, which maximizes return on the power 
spectrum measurement. This would correspond to a detec- 
tion threshold for which the false-positive rate is 50%. In 
this sense, we can choose n — 1 as our selection threshold. 
The selection threshold problem can be dealt with in a more 
sophisticated way. The survey measures a three-dimensional 
map of the sky. Each pixel in that map will have some signif- 
icance of detection. Clearly the large number of low signifi- 
cance pixels do collectively contai n information, i f they can 
be averaged in a meaningful way. Zh ang fc Penl i2005l) de- 



scribe one such algorithm to extract the luminosity function 
deep into the noise. 

With these large samples of galaxies, the galaxy cluster- 
ing can be precisely measured. This all ows constraining cos- 
molo gy through the baryon oscillation jAbdalla fc Rawlingsl 
2005). The redshift distortion of galaxies allows the measure- 
ment of galaxy velocity power spectrum to high accuracy 8 . 
In this paper we do not attempt to utilize all information 
of these galaxies. Our primary goal is to demonstrate the 
feasibility of using future 21cm radio survey to do precision 
lensing measurement. We notice that the properties of lens- 
ing magnification depends strongly on the mass threshold 
or equivalently n. We will explore different n ^ 1 to find 
suitable choice for different quantities. For n > 3, since the 
number of false peaks is much smaller than the number of 
detected galaxies, one can safely neglect all errors caused by 
false detections. But for n < 3, one has to take them into 
account. 



3 CROSS CORRELATIONS OF DIFFERENT 
REDSHIFT BINS 

3.1 Cosmic magnification preliminary 

Cosmic magnification changes the galaxy number overden- 
sity S g to 9 



5„+2((x-1)k + 0(k 



(7) 



Here, a = -/ (> F c )F c /f(F c ), /(> F) is the number of 
galaxies brighter than F, F c is the flux limit adopted and k 
is the lensing convergence. Since (k) = 0, to the accuracy of 
(k 2 ) ~ 10~ 4 , lensing does not change the averaged /(> F). 
Thus a is effectively an observable. 

The usual method to eliminate the intrinsic corre- 
lation of close galaxy pairs is to cross correlate galax- 
ies in two separate redshift bins. Galaxy peculiar veloc- 
ity shifts the position of galaxies in the redshift space by 
cAz < 10 3 km/s. Galaxy correlation becomes negligible 
at scales r > 100/t -1 Mpc, which corresponds to Az ~ 
0.03[ J ff(z)/ J f/ ][r/100/i- 1 Mpc]. Choosing Az > 0.05, resid- 
ual correlations caused by intrinsic galaxy correlation can 
be safely neglected. With accurate measurement of 21cm 
emitting galaxy redshift, this can be done straightforwardly. 

The distribution of foreground galaxies traces lenses of 
background galaxies which cause the magnification effect. 
Thus there exists a correlation between the background mag- 
nification and foreground galaxies. This galaxy-galaxy cor- 
relation is the correlation generally considered in the lit- 
erature. There exists another correlation. Both foreground 
galaxies and background galaxies are lensed by interven- 
ing matter. Thus there exists the background magnification- 
foreground magnification correlation. The combined corre- 
lation is 



{6^(9f,z f )5^(9 b ,z b )) = 2(a b -l){K b 5(8 f ,z f ) 

+ 4( a/ - l)(a b - l){K fKb ) 



(8) 



8 Zhang& Pen, 2005, in preparation 

9 In this expression, we neglect all high order ter ms throu ghout 
this pa per. These terms increase the lensing signal jMenard et alj 
2003a) and thus improve the correlation measurement. 
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Here the subscript / and b denote foreground and back- 
ground respectively. The first term in the right side of the 
equation is the magnification-galaxy correlation and the sec- 
ond term is the magnification-magnification correlation. 
The observed galaxy surface density is 



S 9 = 



n g (z){l 



5 L )dz 



(9) 



Its correlation function is Eq. |H1 weighted by the differen- 
tial galaxy number distribution n g and clustering signal. 
By Limber's approximation, the corresponding 2D angular 
power spectrum of the magnification-galaxy cross correla- 
tion is given by 

(10) 



2tt 



3fi r 



2c- 



± 2 g m ( — ,Zf)G b (zf)n g (zf)xfdzf 
Xf 



Here, Q m is the present day matter density with respect 
to the cosmological critical density. Nf = J n g {zf)dzf is 
the total number of foreground galaxies. A gm = b g r g A 2 m 
is the galaxy-matter cross correlation power spectrum. We 
assume b g r g = 1 (b g is the galaxy bias and r g is the cross 
correlation coefficient between galaxies and dark matter). 
The matter power spe ctrum A^ is calcul ated using the 
BBKS transfer function jBardeen et alll986f) and its nonlin- 
ear evolution is calculate d by Peacock-Dodds fitting formula 
]Peacock fc Dodd'3 [l996') . G b , the kernel of the background 
magnification effect, is given by 



G h (z) 



1 + 2 

N b 



w(x,Xb)n g (z b )2(a(z b ) - l)dz b . (II) 



Here, w(x, Xs) is the lensing geometry function. For the flat 
universe we adopt, the geometry function is simplified to 
w {XiXs) = x(l — x/Xs)- The strength of the magnification 
effect relies on both the strength of lensing, which prefers 
background galaxies with higher z, and a — 1, which prefers 
deeper slope of the mass function at the mass threshold. 

The power spectrum of the corresponding 
magnification-magnification cross correlation is: 



2tv 



\ 2c 2 



7T 

T 



(12) 



Ai n (-,z)G b (z)G f (z) X d X - 
X 

Here, G/, the kernel of foreground magnification effect, is 
given by 

G f (z) = 1±1 ^ f '™* w(x,Xf)n g (z f )2(a(z f ) - l)dz f (13) 



/ 



where Nf is the total number of foreground galaxies. The 
amplitude of relies on both the lensing signal of fore- 
ground and background galaxies. The higher the Zf and z b , 
the stronger the correlation signal. It also depends strongly 
on ctf — 1 and a b — 1. If the mass threshold is larger, a — 1 
is generally larger. 

Depending on the choice of foreground bins and fore- 
ground galaxy selection criteria (which determines a and 
thus the strength of magnification effect), either C^ 3 or 
can dominate. In H3.2I we discuss cases where C MS dominates 
and in ^3.31 we discuss cases where C w dominates. 





Figure 3. The predicted accuracy of cosmic magnification-galaxy 
cross correlation power spectrum C m measured by CLAR 5 year 
survey. Foreground galaxies (0.5 < Zf < 1.0) are selected with 
a 2<r selection threshold. The optimal survey coverage is around 
several hundred square degrees. For such configuration, C M9 can 
be measured to ~ 20% accuracy at I around several thousand 
(lower panel, bin size Al = 0.2/). C MM is several percent of C M9 
and does not show up in the plot. In top panel, we further split 
background galaxies into several redshift bins. This lensing to- 
mography allows the measure of the evolution of matter distribu- 
tion. It also allows the measure of relative change of the comoving 
distance as a function of z to better than 10% in three redshift 
bins at z b > 1. 



3.2 Magnification-galaxy power spectrum 

For a sufficiently wide foreground galaxy redshift distribu- 
tion with low median redshift, the lensing effect and the 
magnification (a — 1) effect are both small. Generally, in 
this case, is much smaller than C^ 9 . In this section, we 
fix the foreground galaxy distribution (0 < Zf < 1.0) and 
vary the background galaxy redshift distribution. For this 
choice of foreground galaxy distribution, C MM is ~ 1% of 
C M9 (Fig. |3| - The correlation signal peaks at I ~ 10 4 , where 
the fluctuation is ~ 10%. 

C^ g is the projection of Ag m along the line of sight (Eq. 
11011 . Given a cosmology, A^ m (k,z) can be extracted using 
the inversion m ethods applied to galaxy surveys and len s- 
ing surveys Ce.g. lDodelson et al] <2002l) : |Pen et al] ]2003rJ) ). 
Ag m (fc, z) contains valuable information of galaxy clustering 
and can put strong constraint on halo occupation distribu- 
tion. 

Cosmological information is also carried in the geometry 
term of lensing, in our case, x m C M9 and C' 1 . By fixing the 
foreground galaxy distribution and varying the background 
galaxy distribution, one can isolate Xb from Ag m and mea- 
sure the dependence of Xb on z b . This method al lows an in- 
dependent and robust c onstraint on cosmology iJain et al] 
120031: IZhang et al]l2003j) . 
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Figure 4. Similar as Fig. |3 but for SKA 5 year survey. Fore- 
ground galaxies are chosen to be at 0.7 < Zr < 1.0. To measure 
C M9 , optimal sky coverage should be around ten thousand de- 
grees. Under such configuration, C m can be measured to several 
percent accuracy at 1000 < I < 10 4 . The dashed line in the lower 
panel is C^/C^ 9 for 6400 deg 2 sky coverage. The dependence of 
comoving distance x{z) on redshift can be measured to be better 
than 1% in 7 redshift bins at Z > 1. 



The statistical error in the C' 9 measurement is 



AC™ 



(C b a + C b shot )(Cl + Ci ot ) 



{21 + 1)AZ/ Sky 



(14) 



Here, Cg is the auto correlation power spectrum of the back- 
ground galaxies, which includes contributions from (SgSg), 
(«(,«(,} and (SgKb). For high redshift background galaxies, 
(KbKb) is not negligible comparing to the intrinsic galaxy 
correlation {8g5g} and thus has to be taken into account. 
Cshot = 47r/ sky [l + N* ois JN b ]/N b is the background shot 
noise power spectrum. The extra factor 1 + N^ oisc /Nb ac- 
counts for the contaminations of false peaks in the selected 
background sample. These false peaks does not correlate 
with each other, so their only effect is to increase the shot 
noise. For selection thresholds above 2-a, false positive rate 
is small and this extra factor can be neglected. Cf is the 
auto correlation power spectrum of the foreground galaxies. 

Statistical errors strongly depend on / sky (when fixing 
total observation time). This dependence is complicated. (1) 
/ak y affects a(zf). Shorter integration time per unit area is 
required in order to survey for a larger sky area. This in- 



10 At frequencies below 1.4 Ghz, radio sources have smooth con- 
tinuu m spectra, so they can be subtracted away in the frequency 
space iWang et al 12005ft . Their only effect is to contribute to T S y S 
and cause a fluctuation in T sys across the sky. This can introduce 
a false correlation, but it is routine in observation to eliminate 
this effect and thus we do not consider it in this paper. 



creases system noise per beam and thus increases the se- 
lection mass threshold of high z galaxies. Since the mass 
function is steeper at higher mass, a(zf) increases and C^ 9 
increases. (2) / s k y affects the relative distribution of galax- 
ies. Larger / s k y survey detects relatively more low z galaxies. 
Since Cf 9 is proportional to A 2 gm weighted by the distribu- 
tion of foreground galaxies and matter clustering is stronger 
at lower z, larger sky coverage tends to increase C^ 9 . But 
on the other hand, the lensing effect is smaller for lower z 
galaxies. This has the effect to decrease w(%, X&) an d thus 
decrease C M9 . Furthermore, the noise terms Cg and Cg in- 
crease. (3) / s k y affects the cosmic variance. (4) / s k y affects 
the total number of foreground and background galaxies and 
thus changes the shot noise. The lower panels of fig. OH & HI 
show the dependence of AC M9 /AC M9 on sky coverage. If 
the sky coverage is too small, the cosmic variance is large. 
If the sky coverage is too large, too few background galaxies 
can be detected. Shot noise begins to dominate at relatively 
large scales. The choice of selection threshold n (n-a) also 
affects the statistical errors. Larger n increases a and thus 
increases C M9 . But it also decreases the number of detected 
galaxies and thus increases shot noise. 

Both CLAR and SKA can measure C M9 precisely (fig. 
For CLAR, the optimal sky coverage is around sev- 
eral hundred square degree. C^ 9 can be measured to ~ 20% 
accuracy for bin size Al = 0.21. It can also measure C M9 at 
several redshift bins and allows isolating geometry. For SKA, 
the optimal sky coverage is around ten thousand square de- 
gree. C^ 9 can be measured to ~ % accuracy. The size of 
background redshift bins can be as narrow as Az < 0.1. 
The change of Xb( z ) can be precisely measured at z > 1. 

The best result from optical galaxy surveys by far is 
an 8a detect ion by SPSS, using a n optimal quasar weight- 
ing function dScranton et aljl2005T) . CLAR will detect much 
less galaxies and cover much smaller sky area, but even 
CLAR can reach ~ 5-a at more than 10 independent bins 
of width Al = 0.21. If adopted similar background galaxy 
weighting function, the result can be further improved. So 
CLAR and SKA can do much better than optical galaxy 
surveys. This is well explained by Eq. 1141 (1) The number 
of z > 1 CLAR galaxies is of the order 10 6 and is much 
higher than the number of SDSS quasars (~ 2 x 10 4 ). This 
significantly reduces the shot noise term Cg hot in Eq. 1141 
(2) CLAR foreground galaxies have broad redshift distri- 
bution, which roughly matches the lensing kernel of back- 
ground galaxies and thus amplifies the cross correlation sig- 
nal. For SDSS, quasars lie at z > 1 while most galaxies 
locate at z <C 0.5. Since these quasars are mainly lensed by 
matter at z > 0.5, the cross correlation signal is weak and 
This is likely the dominant re ason why the 
measured cross correlation s by 2dF and SDSS llMvers et alJ 
12003!: IScranton et aljEooih do not have as good S/N as one 
would naively expect from their large galaxy and quasar 
samples. 

It is difficult for optical surveys to detect large number 
of high z galaxies. But 21cm radio surveys can. This is an 
inherent advantage of 21cm radio surveys to measure cos- 
mic magnification. Furthermore, since 21cm galaxies have 
precise redshift measurement, one can weight foreground 
galaxies deliberatively to optimize the cross correlation mea- 
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Figure 5. The predicted accuracy of C ; MM , magnification- 
Magnification cross correlation power spectrum in two redshift 
bins, by SKA 5 year surveys. The measurement of C ; MM allows 
constraining cosmology and studying matter clustering without 
the complexity of galaxy bias. We choose foreground galaxies at 
1.6 SC Zf ^ 2.0 and background galaxies at 2.0 ^ 25 ^ 4.0. 
We treat Cj 19 as contaminations. Since C MM depends on a — 1 
of foreground galaxies while C M9 does not, we vary the selec- 
tion threshold of foreground galaxies to change C^/C 9 . We 
fix the selection threshold of background galaxies as 4o\ For 6cr 
and lOcr selection threshold of foreground, galaxies, systematic er- 
rors (C M9 ) is comparable to statistical errors (error-bars of data 
points). 

surement. This is another inherent advantage of 21cm radio 
surveys. 

3.3 Magnification-Magnification power spectrum 

The main strength of the cosmic shear power spectrum and 
bispectrum to constrain cosmology lies in the fact that the 
prediction of these quantities only relies on the matter clus- 
tering whose theoretical understanding is robust. The pre- 
diction of C" MM is as straightforward as the prediction of the 
cosmic shear power spectrum, which does not require the 
knowledge of complicated galaxy bias, as C M9 does. So, the 
measurement of C MM would allow robust constraints on cos- 
mological parameters and dark matter clustering, as cosmic 
shear power spectrum does. In this section, we will show 
that SKA is straightforward to measure the magnification 
(foreground)-magnification (background) cross correlation 
power spectrum C MM . 

For the purpose of constraining cosmology using C MM , 
C M9 should be treated as contamination and marginalized 
over. C ; MM depends on the magnification strength of fore- 
ground galaxies, while C^ 9 does not. By increasing the red- 
shifts of foreground galaxies, the lensing signal increases and 
a — 1 also increases, due to higher mass selection threshold 




Figure 6. Auto correlation angular power spectrum of galaxies. 
We disregard close pairs within redshift separation Az < 0.1 and 
thus eliminate intrinsic galaxy clustering. We try different galaxy 
distribution. Lines with error bars, from bottom to top, corre- 
spond to z > 1.5, z > 2.0, and 2 > 2.5, respectively. For our 
choice of 4<r selection threshold and z > 1.5, C MM dominates over 
O* 9 . For higher redshift, the luminosity function is steeper at the 
limiting flux. Larger a — 1 then increases C' 1 with respect to 
C M9 . On the other hand, high redshift galaxies are mainly lensed 
by low z matter distribution. The higher the redshift, the less 
likely that galaxies can be lensed by matter distribution in the 
same redshift bins. This also increase O" 2 with respect to C 9 . 
We assume no evolution in the HI mass function. Realistic evolu- 
tion scenarios would result in more galaxies and thus allow better 
measurement. 

at higher redshift, thus C ; MM increases with respect to Cj 13 . 
Since C 9 is proportional to the strength of matter clus- 
tering, it decreases when increasing redshifts of foreground 
galaxies. But these requirements to increase C MM with re- 
spect to C M9 can be at odds with the requirement to reduce 
statistical errors. For example, increasing the mass selection 
threshold or redshifts of foreground galaxies reduces Nf and 
thus increases shot noise. 

We try different foreground redshift bins, selection 
threshold and sky coverage to optimize the measurement 
of C MM such that C M9 /C MM is smaller or comparable to the 
statistical error. For SKA, it is indeed possible to measure 
C MM and control both systematic errors (C^ 9 ) and statistical 
errors to 10% level (Fig. 0. Since such measurement requires 
large number of high z galaxies, it is extremely difficult for 
optical surveys to realize it. 



4 AUTO CORRELATION 

With the precision measurement of galaxy redshift and large 
amount of high z galaxies, one can extract the cosmic mag- 
nification from the galaxy auto correlation measurement. In 
the correlation estimator, we throw away pairs with redshift 
separation \z\ — Z2I ^ Az c . We choose Az c = 0.1, which cor- 
responds to comoving separation r c ~ 180/i _1 Mpc at z = 1 
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and r c ~ 100/i~ 1 Mpc at z = 2. In Fourier space, this cor- 
responds to cut off the power at k c < l/r c < 0.01/i/Mpc. 
Applying the Limber's equation in Fourier space, the angu- 
lar fluctuation at multipole I is contributed by the spacial 
fluctuation at k — Then an effective cut off k c in Fourier 
space corresponds to an effective cutoff at l c = k c \ < 20. 
So, under the Limber's approximation, one can neglect the 
residual intrinsic correlation of galaxies at Z > 20. 

One can further quantify the residual intrinsic cluster- 
ing. The angular correlation it produces is 



n g {z) 



dz 
dx/dz 



(15) 



x / UVx 2 2 + (A X ) 2 ,z)dA X 



Here, £ 9 is the galaxy correlation function. Numerical calcu- 
lation shows that |wf c (#)| is smaller than 10 _J at all scales. 
For example, |wj C (~ 1 )| is less than ~ and |tOic(~ 1°)| 
is less than several x 10 -7 . One can further convert Wi C (9) to 
the corresponding C%. We find that \l 2 Ci / (2tv)\ < 10" 5 at all 
scales. Specifically, at I ~ 100, \1 2 Ci/(2ty)\ < several x 10~ 6 
and at I ~ 1000, \1 2 Ci/(2ty)\ < several x 10" 6 . The angular 
fluctuation caused by residual galaxy intrinsic clustering is 
roughly 1% of at I ~ 100 and much less than 1% at 
smaller scales. So, the close pair removal procedure effec- 
tively eliminates all intrinsic galaxy clustering. 

The auto correlation function is composed of two parts, 
the one arising from the auto correlation of cosmic mag- 
nification and the one from the cross correlation between 
cosmic magnification and S g . The magnification auto corre- 
lation power spectrum is 



/ 2 Cf 
2tt 



3^l m B 2 
2c 2 



A 2 m (i,z)G(z) 2 fl(z)xdx(16) 
X 



where N is the total number of galaxies in the corresponding 
redshift bin. This expression differs from Eq. 1121 only by a 
factor /| («), which arises from the close pair removal. /2(a) 
is given by 

fl(z) = [^f|] 2 y^(x,Xi)(ai-l)^(x,X2)(a 2 -l) 
x n g (z 1 )n g (z2)0(\z 1 — 22 1 — Az c )dzidz2 . (17) 

The function 9(x) = if x < and 6(a:) = 1 if x > 0. G{z) 
is defined analogous to G b and Gf. 

The power spectrum of the magnification-galaxy cross 
correlation function is 

Z 2 C™ 



2tt 



2c 2 



(18) 



x A 2 gm {— l r -,z)G(z)f 1 {z)n g (z)x{z)dz . 
X( z ) 

The effect of close pair removal is carried by fi (z) 
l + z 



AW = ^ (19) 

X n g (z b )(a(z b ) — l)Q(\z b — z\ — Az c )dz b . 

Again, as explained in 13. 31 for the purpose of constrain- 
ing cosmology, we treat C M9 as contaminations of C M/i . For 
the SKA, both the statistical errors of C MM and systematics 
errors (C' lg ) can be controlled to better than 10% level, if we 



only use z > 1.0 galaxies (Fig.|SJ. For CLAR, the detectabil- 
ity of C MM is sensitive to the evolution of HI mass function. 
Assuming the conservative no evolution model, C MAl can be 
detected at the several a level(fig.|HJl. Though this close pair 
removal method is successful for 21cm emitting galaxies, it is 
essentially infeasible for optical galaxies because of (1) insuf- 
ficient number of high z optical galaxies and (2) inaccurate 
photo-z redshifts and too expensive spectroscopic redshifts. 



5 COSMIC MAGNIFICATION BISPECTRUM 

Lensing bispectrum contains valuable and often complimen- 
tary information on cosmology and the large scale struc- 
ture, comparing to the 2-point correlation power spec- 
trum jBernardeau et alJll997t lHujll999t iBernardeau et al 
£3bJ;[ 



al. 

in 



| 2002tlMenard et alJl2003bl : iPen et al]l2003at iTakada fc Jai: 
2004). But current data only allows low significance 
detection of the skew n ess at several a n gular scales 
iBernardeau et alJ 12001 IPen et al] l2003al : iJarvis et all 
l200^) ~21cm radio surveys can do much better. In this sec- 
tion, we will show its feasibility by focusing on the bispec- 
trum of galaxies in the same redshift bins. We throw away 
close pairs with \zi — Zj\ < 0.1 where Zi (i = 1,2,3) is the 
redshift of each galaxies. 

The bispectrum comes from four parts, {fififi), (fifig), 
(ngg) and (ggg). The (figg) and (ggg) terms are negligi- 
ble following similar argument in J21 The (nflfi) term con- 
tributes a bispectrum 



3f2 m -Ho 
2c 2 



G 3 (z)f!(z) X - 



xBs( — , — , — ;x)dx 
XXX 



(20) 



Here, fs(z) takes the effect of close pair removal into 
account. Bs(k\,k2,k?,\x) ls t ne matter density bispec- 
trum, which is calculated adopt ing the fitting formula of 
IScoccimarro fc Couchmanl i200lTl . 

The (fifig) term contributes another bispectrum 



/ 3£l m H 2 
\ 2c 2 



G 2 {z)f i {z) x - i 



, n f h h h ^n g (z) 



(21) 



Here, }a{z) takes the effect of close pair removal into ac- 
count. We explicitly show the galaxy bias b g in the above 
equation, though we adopt b g = 1 in the estimation. The 
factor 3 comes from the permutation of {fj,ifijgi„}. 
The shot noise of the bispectrum is 



-Bshot^l, I2, I3) — C N + Cn{Ci + C2 + C3) 



(22) 



Here, Cn = 47r/ s k y /A i ' is the shot noise power spectrum. 
Full evaluation of the sample variance of bispectrum involves 
integrating 6-point nonlinear density correlation function. 
Since we do not have robust theoretical predictions or simu- 
lation results of 6-pt nonlinear correlation function, we only 
consider the Gaussian sample variance. The statistical error 
of corresponding bispectrum .B123 is 



ABi 



(C 1 C2C 3 + B 2 hot ) 



Nr 



(23) 
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Figure 7. Forecasted measurement of bispectrum by SKA. We 
assume no evolution in the HI mass function. The solid lines are 
^fiMM °f equilateral configuration (£1 = I2 = I3), the dot lines are 
B/j,fx g and the dash lines are the statistical error. For statistical 
error, we adopt the I bin size Al = 0.21 and the angular bin 
size 10°. We try different galaxy redshift ranges and selection 
thresholds. The top left, top right, bottom left and bottom right 
panels are the results of (3a, z > 1.5), (5cr, z > 1.5), (3<x, z > 1.0) 
and (5<t, z > 1.0), respectively. 



Here, C» = C^{U) + C M {U). N123 is the number of inde- 
pendent combination of h,h,h used to obtain the averaged 
B123. For a rectangle survey area with x axis size Q x and y 
axis size 9 y , independent modes are 1 = (2nm/8 x ,2im/9 y ), 
where m,n — 0, ±1,±2, There is a constraint that 
Z3 = —li — I2, then the total number of independent combi- 
nation is 

^123 = dmi x dni x drri2xdn2y (24) 

= (j^7^) hdhhdhl-Kdd^ 

where #12 is the angle between li and I2. We only show the 
result of li ~ I2 — I3, for which, we choose Al = 0.2Z and 
A0i2 = tt/18 (10°). 

At large scales (I < 1000), cosmic variance prohibits 
the measurement of -B M/J(J and B^^g. But at scales / ~ 10 4 , 
either -B MMM or 5 M/Jg or both can be measured to be better 
than 10% accuracy by SKA (Fig. 0. As expected, higher 
z and/or higher selection threshold result in higher lensing 
signal and amplify B^^ with respect to S MM9 . The com- 
putation of three point function may appear computation- 
ally challenging, requiring the enumeration of A?" 3 triangles 
with N ~ 10 8 . Recently, linear algorithms have been devised 
which resolve these problems JZhang fc Pen! 120051) . Again 
the requirements of precise redshift measurement and high 
z galaxies prohibit the measurement of lensing bispectrum 
through cosmic magnification of optical galaxies. 



Figure 8. C MM and C M9 that would be measured in the same 
redshift bin by CLAR. We try different evolution models, as ex- 
plained in JH] Top left panel assumes no evolution and we choose 
galaxies selected above 2a at z > 1.5. Other panels assume evo- 
lution models explained in the text. Top right panel uses galaxies 
above 3cr at 2 > 1.5. Bottom right panel uses galaxies above 2a 
at 2 > 2.0, Bottom left panel uses galaxies above la atz > 2.5. 



6 EVOLUTION EFFECT 

We have demonstrated the feasibility of measuring lensing 
power spectra and bispectrum in cross correlation of galaxies 
in two redshift bins and in auto correlation of galaxies in the 
same redshift bin. We caution on that these results (fig. QUI 
I3I4I5I & I6I7H should be regarded as conservative estimation of 
the power of radio surveys to measure cosmic magnification. 
There are several reasons. One is that in this paper we only 
use galaxies in certain redshift ranges and above certain se- 
lection threshold. Better measurement of and C MS can 
be obtained by cross correlating galaxies above la detec- 
tion threshold at all redshifts. To estimate how much can 
one gain requires the design of careful weighting on different 
redshifts and luminosity (HI mass). This work is beyond the 
scope of this paper. Another reason is that we have assumed 
no evolution in the HI mass function. Evolution effect is very 
likely to improve the accuracy of lensing measurement by 
providing many more detected galaxies. We will investigate 
the evolution effect in this section. 

As discussed in J|H the observed Qmh shows a factor 
of 5 increase from z = to z ~ 3. Its evolution can be ap- 
proximated as g(z) = (1 + z) 2 ' 9 exp(— z/1.3). Thus we have 
a constraint of no(z)M*(z) = no(z = 0)M*(z — 0)g(z). 
There is little solid constraint on the evolution of no or M* 
separately. But the observation of damped Lyman-a sys- 
tems and Lyman-Limit systems provides some indirect con- 
straints. Damped Lyman-a systems have HI column density 
A^hi ) 2 x 10 20 cm -2 . If the size of the corresponding HI 
regions is ~ 30 kpc/h, then the total HI mass is ~ 10 10 Mq. 
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Thus these damped Lyman-a systems are likely part of cor- 
responding massive HI (proto-)galaxies. On the other hand, 
Lyman-Limit systems have much smaller HI column den- 
sity and are likely part of less massive HI galaxies. The ra- 
tio of damped Lyman-a systems abundance with respect to 
Lyman-limit systems decreases after z = 3. This implies that 
there may be fewer massive HI galaxies after z ~ 3 and thus 
an evolution of M*{z). 

Since the constraint to either no or M* is weak and it 
is likely that both no(z) and M*(z) evolve, we explore three 
evolution scenarios. (A) No evolution in M,(z). no(z) = 
n (z = 0)g(z). (B) No evolution in n (z). M*(z) = M»(z = 
0)g(z). (C). n (z)/n (z = 0) = M*{z)/M*[z = 0) = g{zf/ 2 . 

The number of z > 1 galaxies increases by at least a 
factor of 5 for these evolution scenarios. Taken the evolution 
effect into account, even CLAR can measure C MA1 to ~ 10% 
accuracy (systematic and statistical, fig.|HJ. 



7 DISCUSSION 

We further address that the results shown in this paper only 
utilize a small fraction of cosmic magnification information 
contained in 21cm emitting galaxy distribution. (1) We only 
tried several bins of galaxy redshift distribution and sev- 
eral galaxy selection threshold to demonstrate that cosmic 
magnification can be measured to high accuracy. To utilize 
the full lensing information, one needs to divide galaxies 
into many redshift bins and selection threshold bins. One 
then measures (N-point) auto correlation functions of each 
bins and cross correlation functions between different bins. 
In principle, one can develop optimal weighting scheme to 
combine all measurements to get the best measurement of 
lensing statistics and lensing-galaxy statistics. (2) We did 
not attempt to separate the cosmic magnification auto cor- 
relations and the cosmic magnification-galaxy correlations. 
These two classes of correlation have different dependence 
on the selection threshold of galaxies. These dependences 
are straightforward to predict and can be applied to sepa- 
rate two components. Such component separation improves 
the robustness of constraining cosmology and large scale 
structure significantly. The cosmic magnification auto cor- 
relations an d the geometry of cosmic m agnification-galaxy 
correlations (Jain et all2003l IZhang et all 2003) are ideal to 
constrain cosmology and matter clustering. The amplitude 
and angular scale dependence cosmic magnification-galaxy 
correlations are ideal to constrain halo occupation distri- 
bution. Advanced analysis methods are required to address 
the above two issues and to utilize full information of cosmic 
magnification information in 21cm emitting galaxy distribu- 
tion. (3) We note that individual galaxies can be resolved 
with SKA, this allows the measurement of cosmic shear. It 
may also allow an independent measure of cosmic magnifi- 
cation. SKA resolution allows the measurement of an incli- 
nation angle. If galaxies at high redshift also follow a Tully- 
Fisher relation, the lensing effect can also be large compared 
to shot noise and can be extracted. Since we do not know the 
evolution of the scatter in the Tully-Fisher at high redshift, 
we do not use this information in this paper. It is likely that 
real surveys can do significantly better than our estimates. 

Utilizing all information of cosmic magnification infor- 
mation in 21cm emitting galaxy distribution, the relative er- 



ror of cosmic magnification will be much smaller than what 
shown in this paper. Upon this precision era, one needs to 
improve the theoretical prediction to better than 1% accu- 
racy and at the same time, understand possible 1% system- 
atics. 

Our prediction of cosmic magnification is simplified 
in two ways. (1) We only considered the leading order 
term of cosmic magnification (Eq. 0. Higher order terms 
are known to be cap able of generating 10% effect (e.g. 
iMenard et all J2003ali l and have to be included to inter- 
pret data at the forecast accuracy. (2) For the lensing con- 
vergence, we neglected high order corrections caused by 
lens-lens coupling and deviation from Born's approxima- 
tion. These hig h order corrections are known to have several 
percent effect fechneider et alJ Il99gl: IVale fc White! 120031: 
iDodelson fc Zhandl2005F and should be taken into account. 

Source-lens coupling llBernardeaulll998l : lHamanall200lT) 
can have non-negligible effect to cosmic shear. It arises from 
the fact that measured cosmic shear is always weighted by 
the number of observed galaxies, which also trace the matter 
distribution. But this effect does not exist in cosmic magni- 
fication where we directly correlate the numbers of observed 
galaxies at two different redshifts and directions. 

Several other approximations in our cosmic magnifica- 
tion measurement only introduce negligible corrections. (1) 
we have assumed that the luminosity function /(> F) is 
the same everywhere and thus a is the same everywhere. 
This picture is over simplified. /(> F) can have environ- 
mental dependence. Thus in principle (aa) ^ (ex) 2 - This af- 
fects the prediction of correlations where two or more cosmic 
magnification terms present (e.g. in C MM , -B MMM and -B M(JS ). 
But this effect is very small, since a only depends on lo- 
cal environment, for two redshift bins with modest sepa- 
ration Az > 0.05, (0102) — (qi)(q2)- The close pair re- 
moval procedure further guarantees that even for the same 
redshift bin, (0102) — (ai)(a2). Thus, one can safely ne- 
glect this effect. (2) Residual intrinsic clustering causes 

~ 1% correction at I ~ 100 O • Since cosmic variance 
1/2 

at I ~ 100 is > 0.01/ sk y' . One needs to worry about this 
effect only for full sky surveys. Furthermore, it can be re- 
duced by extrapolating galaxy correlation function £ s at 
smaller scales measured from the same survey to relevant 
scales (> 1007i _1 Mpc). 



8 CONCLUSIONS 

We have made simple forecasts for future radio surveys to 
measure gravitational lensing. We found that radio surveys 
can be precise sources for lensing measurements, and that 
lensing magnification is measurable because redshifts are 
known and many galaxies can be detected. CLAR and SKA 
are expected to measure the dark matter power spectrum 
and galaxy-matter cross correlation to high accuracy. The 
estimates are conservative, since evolution in HI mass func- 
tion, as suggested by observations, will improve the above 
results significantly. Also, many effects will increase the sen- 
sitivity. Unfortunately, these effects, which include Tully- 



1 But it causes -C 1% correction at smaller scales. 
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Fisher relations, are difficult to quantify at high redshifts, 
for which we neglect in this paper. 

More complete statistical information is available at 
lower signal to noise levels and using non-Gaussian statis- 
tics. We have made estimates of the three point statistics, 
which appear promising, and are expected to improve the 
information that can be gained by lensing. 
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